Correlation energy of the spin-polarized uniform electron gas at high density 
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The correlation energy per electron in the high-density uniform electron gas can be written as 
Ec{ra,C) = Ao(C)liir's + eo(C) + \i{C,)rs\nrs + 0{rs), where rs is the Seitz radius and C, is the 
relative spin polarization. We derive an expression for Ai((") which is exact for any C,, including the 
paramagnetic and ferromagnetic hmits, Ai(0) and Ai(l), and discover that the previously published 
Ai(l) value is incorrect. We trace this error to an integration and limit that do not commute. The 
spin-resolution of Ai((") into contributions of electron pairs is also derived. 
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The final decades of the twentieth century witnessed 
a major revolution in solid-state and molecular physics, 
as the introduction of sophisticated exchange-correlation 
models^ propelled density-functional theory (DFT) from 
qualitative to quantitative usefulness. In principle, the 
foundation of DFT is the Hohenberg-Kohn theorerrP but, 
in practice, it is usually the supposed similarity between 
the electron density in a real system and the electron 
density in the hypothetical uniform electron gas (UEG).'^ 

The three-dimensional UEG is characterized by a den- 
sity p = P'l -\~ Pi, where p^ and p^ is the (uniform) density 
of the spin-up and spin-down electrons, respectively. In or- 
der to guarantee its stability, the electrons are assumed to 
be embedded in a uniform background of positive charge. 

In 1965, Kohn and Sharn^ showed that the knowledge 
of a analytical parametrization of the UEG correlation 
energy allows one to perform approximate calculations for 
atoms, molecules and solids. This led to the development 
of various spin-density correlation functionals (VWN,'^ 
VTif^ PW92,'2j etc.), all of which require information on 
the high- and low-density regimes of the spin-polarized 
UEG, and are parametrized using results from near-exact 
Quantum Monte Carlo (QMC) calculations .l^ti^ 

However, inspired by Wigner's seminal work,'^^ Sun, 
Perdew and Seidl have recently shown that the correlation 
energy of the UEG can be estimated accurately without 
any QMC input They used a "density-parameter in- 
terpolation" (DPI) between the (near-)exact high- and 
low-density regimes, which precisely reproduces the first 
few coefficients of the high- and low-density energy 
expansions.!^ Knowledge of these coefficients, of course, 
is essential for such interpolations and is the motivation 
for the present work. 

The high-density expansion of the correlat ion energy 
per electron (or reduced energ y) of the UEG ipHlIIlHMl 

E,{r,,0 = Ao(C) lnr,+eo(C)+Ai(C) In r3+0(r,), (1) 

where Ts = (47rp/3) is the so-called Seitz radius, and 
C — (Pt ~ Pi)/ P the relative spin polarization. It is 
clear that Ao(C), £o(C): '^ilOi • ■ • must be even functions. 
We use atomic units throughout. 

The coefficient Ao(C) can be obtained by the Gell-Mann- 
Brueckner resummation technique,'^ which sums the most 
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FIG. 1. The five spin-scalings as functions C,. 

divergent terms of the series ([T]) to obtain 
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and kf^i — {1 ± (^) is the Fermi momentum of the spin- 
up or spin-down electrons. The paramagnetic^ (C = 01 
and ferromagneticPH (^ = 1) limits are given in Table I 
and the spin-scaling function 
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was obtained by Wang and Perdew.l^ 

The coefficient eo(C) is usually written as the sum 

soiO^eUO+s"^ (6) 

of a RPA (random-phase approximation) term £q{C) and 
a first-order exchange term Eq. The RPA term £q(C) is not 
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TABLE I. Energy coefficients and spin-scaling functions of tlie paramagnetic = 0) and ferromagnetic (C = 1) states of the 
liigli-density UEG. Note tliat a — (97r/4)~^''^ and z{ri) is the Riemann zeta functionP^ 



Term 



Coefficient 



Paramagnetic limit 
£(0), A(0) 



Ferromagnetic limit 
£(1), 



Spin-scaling function 
T(C), A(C) 
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Eq. ([5| 
Ref.lM] 
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TABLE II. Reduced correlation energy ~Ec{ra,l) for the 
ferromagnetic state of the UEG for various rs. 



where 
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0.0154(1) 
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0.0152 


10 


0.0105(1) 




0.0102 




0.0103 


20 


0.00678(2) 


0.00663 


0.00664 


50 


0.00355(1) 


0.00350 


0.00350 


100 


0.002073(3) 


0.002055 


0.002055 



^ Benchmark QMC results taken from Ref. 8 The digits in 

parenthesis represent the error bar in the last decimal place. 
^ Results taken from Ref. 1171 using the DPI (density-parameter 

interpolation) formula with Ai(l) = 0.003125. 

Results from the present work using the DPI formula with 

Ai(l) = 0.004792. 



known in closed form, but it can be computed numerically 
with high precision.'2Sl Its paramagnetic and ferromagnetic 
hmits are given in Table |l] and the spin-scaling function 
Tg(C) = eg(C)/eo(0) can be found using Eq. (20) in Ref.EHl 
The first-order exchange terrrPU jg given in Tableland, 
because it is independent of the spin-polarization, the 
spin-scaling function Tg(C) = eo(C)/£o(0) = 1 is trivial. 
The coefficient Ai(C) can be written similarlji^i as 
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are the RPA and second-order exchange contri buti ons 
and a = {9tt/4)~^^^ . The integrand functions ar^^E^ 
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Ai(C) = A?(C) + A?(C), 
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Carr and Maradudin gave an estimat^^l of Ai (0 ) and this 
was later refined by Perdew and coworkers .l^EIl 

However, we have found that the integrals in Egs. [8]and 
[9] can be evaluated exactly by computer software,'^" giving 
the paramagnetic and ferromagnetic values in Table [T] and 
the spin-scaling functions 
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where Li2 is the dilogarithm function!^ 

The spin-scahngs Ao(C), Tg(C), T^{C), Af(C) and A^rn 
are shown in Fig.[l] highhghting the Hoffmann minimunPSl 
in Tg(C) near C = 0.9956 and reveaiing a simiiar minimum 
in Af(C) near = 0.9960. Such minima seem to be 
ubiquitous in RPA coefficients. 

The data in Table |T] yieid the exact values 
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and it is revealing to compare these with recent numerical 
calculations. The estimate Ai(0) w 0.0092292 by Sun et 
al^ agrees perfectly with Eq. ([Is]) but their estimate 
Ai(l) « 0.003125 is strikingly different from Eq. ([l9|. 
How can this discrepancy be explained? 

Following Gell-Mann and Bruecknerp^l and Uedap^ 
Misawa arguect^Sj that the C = and C = 1 limits of the 
RPA and exchange contributions to the correlation energy 
are related by 

E:{r^,l) = lE^,{2-^^\s,0), (20) 
E^{r,,l)^E^{2-^/\,,0), (21) 

and, from these relations, Perdew and Wang inferred^ 



A^^(l) = 2-^/3 A^^(O), 
A5'(1) = 2-4/3 A!j'(0). 



(22) 
(23) 



These are also obtained if the C 1 limit of the integrands 
in Eqs.[8]and|9]is taken before integrating over u. 

Numerical evaluations of Eq. ( [9| a nd analytical results 
from Eq. ([Tt]) confirm that Eq. ( |23[ ) is correct. However, 
numerical evaluations of Eq. [8]and analytical results from 
Eq. (161 agree that Eq. ([22]) is wrong and that, in fact. 
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The error in Eq. ( 22 ) arises from the non-commutivity of 
the C ~^ 1 limit and the u integration, which is due to the 
non-uniform convergence of TZ\{u, Q. 



To show this particular point, let us define 
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It can be easily shown that it is not possible to find 
a function D{u), which is integrable with respect to u 
and dominates Anf{u,C), i-e. V(u,C), \Anf{u,C)\ < 
D{u). Thus, according to the dominated convergence 
theorem, one cannot show that limit and integration can 
be interchanged, and 



lim A7e?(u, C) = 0, 



AA'^(l) = 0. (27) 



However, substituting t = u/k^ in Eqs. (25) and (26), 
one immediately finds a function D(t), which is integrable 
with respect to t and dominates A7?.f(i, C)- This ensures 
the possibility of interchanging limit and integration. It 
yields 
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dt 



(28) 
(29) 
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= 0.00166727, 



which is exactly the difference between the two values of 

A!(i). 

The effect of the coefficient Ai(l) on the reduced corre- 
lation energy of the ferromagnetic state has been studied 
by varying its value in the DPI formula proposed by Sun 
et al. in Ref. [T71 The results have been compared with 
the benchmark QMC calculations of Ceperley and Alder.^^ 
As shown in Table [TH the new value of Ai(l) derived in 
the present study systematically improves the accuracy 
of the DPI correlation energy, especially for small r^. 

In some cases, it is of interest to resolve Ai(C) into 
contributions due to tti and ti electron pairs, such as 

Ai(C) = >^fHO + Ai'^-^(C) + Ai'^-^(C), (30) 
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where i = a or b, and aa' — Y[, or t-i-. Using ( 16 1 and 
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(17), we find 
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The remaining contributions can be obtained using the 
relations 



A 



, (C) = A1^"(-C), 

Ai-^^(c) = i-Ai-^^(c)-Ai-^^(c), 



(34) 
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and are represented in Fig. [2] 

In conclusion, we have found a closed-form expression 
for the coefficient Ai(C) of the r^lnrs term in Eq. ([I]). 
It is valid for any value of C and, in particular, for the 
paramagnetic {( — 0) and ferromagnetic {( — 1) limits. 



This reveals that an earlier derivation of the ferromag- 
netic limit Ai(l) was incorrect because of an inadmissible 
interchange of a limit and an integral. The present result 
has no direct impact on the quantum phase diagram of 
the UEG, because the effect of the coefficient Ai(C) is 
more pronounced in the high-density limit (0 < Ts < 2), 
where the paramagnetic fluid is significantly more sta- 
ble than the ferromagnetic one.l^ Preliminary results on 
higher-order coefficients reveal that they behave similarly, 
and special care has to be taken in future studies. We 
believe that these new results will be useful in the future 
development of exchange-correlation functionals within 
DPT. 
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